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ABSTRACT 

We considered the quark-meson coupling model for a nucleon. The model 
describes a nucleon as an MIT bag, in which quarks are coupled to the scalar 
and the vector mesons. A set of coupled equations for the quark and the meson 
fields are obtained and are solved in a self-consistent manner. We show that 
the mass of a dressed MIT bag interacting with a- and w-meson fields differs 
considerably from the mass of the free MIT bag. The effects of the density- 
dependent bag constant are investigated. The results of our calculations imply 
that the self-energy of the bag in the quark- meson coupling model is significant 
and needs to be considered in doing the calculations for nuclear matter or finite 
nuclei. 



1 Introduction 

For more than a decade the success of quantum hadrodynamics has been 
rather impressive in describing the bulk properties of nuclear matter as well 
as the properties of finite nuclei. The model is rather simple with only a 
few parameters, and yet it has been successfully applied to a great number of 
problems for nuclear matter and nuclear structure. In this model the relevant 
degrees of freedom are nucleons and mesons, but the nucleons having a com- 
posite structure are treated as Dirac particles. n Several years ago a model to 
remedy this problem was proposed by Guichona. He proposed a quark-meson 
coupling (QMC) model, in which quarks and mesons are explicitly dealt with. 
The model describes nuclear matter as non-overlapping, static, spherical MIT 
bagstrD interacting through the self-consistent exchange of scalar (a) and vector 
(ui) mesons in the mean-field approximation. Using this model, he investigated 
the direct quark degrees of freedom in nuclear matter. The model was refined 
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later to include the nucleon Fermi motion and the center-of-mass corrections 
to the bag energylj and was applied to a variety of problems [6-10]. 

In the MIT bag model Bi the bag constant B and a phenomenological 
parameter Z are fixed such that the nucleon mass of 939 MeV is reproduced 
for some bag radius R. The bag constant B produces the pressure to make a 
bubble in the QCD vacuum. Z is to account for various corrections including 
the zero-point motion. In the QMC model for nuclear matter the quark-cr 
coupling constant and the quark-w coupling constant (g^ ) are additionally 
introduced. Usually, B and Z are fixed first so that the mass of the free MIT 
bag becomes equivalent to the nucleon mass of 939 MeV for a certain bag 
radius. We denote by B^ ree and Z^ ree , respectively, the bag constant B and 
the parameter Z for the free MIT bag. B^ ree and Z^ ree thus obtained are 
used in the QMC model calculations, and then g% and g^ are determined to 
reproduce the binding energy per nucleon (B.E. = —16 MeV) at the saturation 
density (p° N = 0.17 fm~ 3 ). 

However, due to the interaction between the quarks and the a- and the u>- 
mesons the mass of a single dressed MIT bag in free space, when calculated with 
gfree anc j ^free d esc rirjed above, may be different from the nucleon mass of 
939 MeV. Such a possible deviation of the mass of a dressed bag from the value 
of 939 MeV has been neglected in previous QMC model calculations for nuclear 
matter. If the deviation in mass is significant, it is necessary to modify the 
parameters B^ ree and Z^ ree before implementing them in the nuclear matter 
calculations. In this paper we investigate this change in the nucleon mass. In 
Section 2 the QMC model for a nucleon is described. Some numerical results 
are presented in Section 3. Section 4 contains a summary. 



2 Quark-meson coupling model for a nucleon 

The Lagrangian density for the MIT bag in which quark fields ip are coupled 
to the a- and u>- fields may be written as 

£ = \ WW ~ M)-f1>) - $ {gll^ + (m, - g«a)) i/> - B(a) 

where F^ v = d^oj v — D^lu^, 9 V is the step function for confining the quarks 
inside the bag, and A s is the (^-function at the bag surface. Here, we assume 
the density-dependent bag constant introduced in Ref. [10] and express the 
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bag constant B{a) as 



B 
B^ 



A a 



(2) 



with M% = 939 MeV. We will neglect the isospin breaking and take m q = 
(m u + md)/2 hereafter. For actual numerical calculations m q will be taken to 
be zero. The mass of a (u>) is taken as 550 (783) MeV. From this Lagrangian 
density the equations of motion for quark fields ip, sigma fields a, and omega 
fields uj^ follow. For ip we have 

[V (i0„ - gluj,) - (m, - g<a)] i>6 v = i(l - i^n^A s , (3) 

where A s = —n ■ d(9 v ) is used. The left hand side of Eq. (||) gives us the 
equation for the quarks inside the bag, and the right hand side gives us the 
linear boundary condition at the bag surface. The equations for a and lo^ are, 
respectively, 



dB(a) 



da 



and 



(4) 
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If we consider only the ground state quarks in the static spherical MIT bag 
and keep only the time-component of lj^, Eq. (^) for r < R becomes coupled 
linear differential equations for g(r) and f(r); 



#) 
dr 

dg{r) 
dr 



2 l^-L + E-(r)g(r) 



E+(r)f(r), 



(6) 
(7) 



where g(r) and f(r) are the radial parts of the upper and the lower components 
of ijj, i.e., 



qt/R 



9{r) 
-id ■ f/(r) 



Xq 



4ir 



where a is the Pauli spin matrix and \q is the quark spinor. Also, 



(8) 



£ + « = j-^o(r) + K-flMr)) 
E-(r) = | - g*u (r) - (m, - g%a{v)) 



(9) 



with tr(r) and ujo(t) being the a- and the time component of the cj-fields. For 
the ground state Eqs. (JsJ) are the equations in the radial coordinate r only. The 
linear boundary condition from the right hand side of Eq. (^), when rewritten 
by using Eq. (g), becomes 

at the bag surface, r — R, and determines the eigenvalue x of the quarks. e q 
is then given by \/ x 1 + (Rm q ) 2 — x for m q = 0. g{r) and /(r) are normalized 
such that the quark density p q {— V>7o ip) integrated over the bag is the unity. 
In the static spherical approximation Eqs. (Q) and (^) are reduced to 



dB(a) 



e(R-r) (11) 



(V 2 -m 2 )a(r) 
and 

(V 2 - ro*)wo(r) = - 5 «(3p,) - r), (12) 

respectively, where p s (= V'VO is the scalar density and 3 is explicitly multiplied 
by the densities (p s and p q ) here to account for the sum over 3 quarks. Equa- 
tions ( |TT| ) and (jl^), which are nothing but the Klein-Gordon equations with 
source terms, may be readily solved using the Green's function defined by the 
equation 

(V 2 -m?)Gi(r,r') = 6{r - r') (i = <r,w). (13) 
For the s-wave state C?,(r, r') may be written asQ 

GAr, r') = ~ — - sinhirmrA e -" l * r > . (14) 

47r rriirr 

Equations (||), (^), (O), and (12) form a set of coupled equations for cr, and 
ujq, which need to be solved self-consistently. 

By solving these equations we can obtain the eigenvalue of the quarks and 
the energy of the nucleon bag. One way to calculate the energy {Em) of the 
nucleon is computing the energy-momentum tensor and using the relation 



E N = J d 6 r Too. (15) 

Too may be written as 

Too = (f^V + B(a)j 6{R-r)-\ ((Vc ) 2 + (m^ ) 2 ) + i ((Va) 2 + (m CT a) 2 ) , 

(16) 
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where E q is given by 

^ = 3^2 --=3--— (17) 
R R R R R 

with the sum over 3 quarks taken into account. Correcting for spurious center- 
of-mass motion in the bag, the mass of the bag at rest is taken to be til 



Mn = y E N - (pc.m.)i (18) 

where (p^.m.) = Sfe=i(Pfc) = 3(x/R) 2 . By calculating Mm for each bag radius 
and minimizing Mjv with respect to the bag radius R, we can get the nucleon 
mass and the bag radius. 



3 Results 

We first present some results of our calculations of the nucleon mass when 
we take the model parameters from Ref. [10]. There are several sets of pa- 
rameters used in Ref. [10]. Here we only consider the cases when R = 0.6 
fm and 8 = 8. In Fig. 1 Mm is plotted as a function of the bag radius R. 
The solid curve represents the free MIT bag mass M N ree without a coupling 
of the quarks with the mesons; g% = g^ = g^ = 0. B^ 4 and Z are taken 
to be 188.1 MeV and 2.030&E2I, respectively, which produce the minimum of 
M N ree at 939 MeV and R = 0.6 fm. Now we include the couplings between 
the quarks and the mesons. We take the parameter sets used to plot Fig. 1 
of Ref. [10]. When g% = and 8 = 8, we have g% = 8.45 and gl = 3.00 
from Table 1 of Ref. [10]. Mm obtained with these parameters are plotted 
by the dashed curve in Fig. 1. Mm is reduced to 848 MeV with R = 0.694 
fm. When g% = 1.0 and 8 = 8, and thus g% = 7.26 and g% = 3.08 from 
Table 1 of Ref. [10], we get Mm as plotted by the dotted curve in Fig. 1 with 
M N = 834 MeV at R = 0.674 fm. Let us consider the case when g% = 2.0 
and 8 = 8. Then from Table 1 of Ref. [10] we have g% = 5.65 and gl = 2.85. 
Mm obtained with these parameters are plotted by the dash-dotted curve in 
Fig. 1 with M N = 830 MeV and R = 0.659 fm. Fig. 1 of Ref. [10] shows that 
for these three different sets of parameters the nuclear matter binding energies 
calculated by the QMC model are the same particularly near the saturation 
density. However, the present calculations show that the nucleon masses cal- 
culated with these parameter sets are considerably different from each other, 
all being far from M N ree = 939 MeV. The couplings induce large attraction 
and reduce the nucleon mass by about 100 MeV. Fig. 1 also shows that the 
attraction is not uniform with respect to the bag radius. There is more attrac- 
tion at larger bag radii, and as a result the bag radius increases significantly. 
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Fig. 1 : Nucleon mass Mjv plotted as a function of the bag radius R for different coupling 
constants, g%, gl and gj? . b]J A = 188.1 MeV and Z = 2.030 are used for all calculations. 
For the dashed curve, g^ = 0, g^ = 3.00 and = 8.45. For the dotted curve, gl = 1, 
g* = 3.08 and g§ = 7.26. For the dash-dotted curve, g% = 2, gl = 2.85 and = 5.65. 

To see how these results come about we consider some simple cases. Let's 
go back to the free MIT bag case; g q a =gl=g® = 0. We still use B 1 ^ = 188.1 
MeV and Z = 2.030. The free MIT bag mass is plotted again in Fig. 2 by 
the solid curve. Let us introduce the coupling between the quarks and the 
w-meson by use of g% — and = 2 and use a constant bag constant, i.e., 
g^ = 0. Then the dashed curve in Fig. 2 is obtained, giving us Mjv = 957 
MeV at R = 0.595 fm due to the repulsive interaction. When g q a = 2 and 
g% = 0, i.e., when there is only quark-cr coupling, attraction is induced and we 
get M N = 926 MeV at R — 0.601 fm, as shown by the dotted curve. When 
we take g% = g^ = 2, cancellation between the repulsion and the attraction 
takes place, and a small repulsion is left as plotted by the dash-dotted curve, 
resulting in M N = 948 MeV at R = 0.598 fm. 

Now we consider the density-dependency in the bag constant using Eq. 
(||). We start with the case that g% = g% = 2 and g^ — 0, corresponding to 
the dash-dotted curve in Fig. 2. This curve is plotted again in Fig. 3 by the 
dash-dotted curve. (M^ ee is also plotted in Fig. 3 by the solid curve.) To see 
the effects of the density-dependency of the bag constant, we take g^ = 1, 
keeping g% — = 2. The calculated Mn is represented by the dashed 
curve. Reducing the bag constant by introducing the density-dependency 
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Fig. 2 : Nucleon mass Mff plotted as a function of the bag radius R for different coupling 
constants, g%, gZ and g£ . Bf/ 4 = 188. f MeV and Z = 2.030 are used for all calculations. 
For the dashed curve, g^ = and g^ = 2 are used. For the dotted curve, = 2 and gZ = 0. 



For the dash-dotted curve, g% 
used. 



2. For all the calculations shown here g^ = are 
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Fig. 3 : Nucleon mass Mjv plotted as a function of the bag radius with the coupling constant 
given in the figure. B^ 4 and Z are the same as in Fig. 2. 



with = 1 causes attraction. For = 3, more attraction is induced, 
as shown by the dotted curve in Fig. 3. This can be seen from Eqs. (11) 
and ([!(]). The source function for the a field in Eq. (|ll|) increases (nega- 
tively) due to dB(a)/d(T, so the a field increases, which results in enhanced 
attraction in the eigenvalue. In Fig. 4 we show the er(r) field for R = 0.6 
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fm as a function of the radial coordinate r. The dotted, dashed, and dash- 
dotted curves in Fig. 4 are obtained with the corresponding parameters used 
for the dotted, dashed, and dash-dotted curves in Fig. 3, respectively. The 
figure shows that as becomes larger the a field increases due to the in- 
crease of the source function. Hence, the smaller nucleon mass for the larger 
gf. Of course, a smaller B(a) in Eq. (16) causes a smaller mass, too. 
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Fig. 4 : a fields for R = 0.6 fm. The coupling constants used here for the dotted, dashed, 
dash-dotted curves are the same as those for the dotted, dashed, dash-dotted curves, respec- 
tively, in Fig. 3. 



As we mentioned earlier, when g^ is large enough, not only the nucleon 
mass becomes smaller, but also the bag radius becomes larger. When the bag 
constant is kept constant as in Fig. 2, the bag radii are essentially fixed as 
R = 0.6 fm. Even if we change g% and g^ to some extent, the nucleon mass 
curves move more or less uniformly in energy. However, when we introduce 
the density-dependency in the bag constant, the bag radius increases as shown 
in Figs. 1 and 3. The density-dependent bag constant brings the Mjv curve 
downward in energy, but not uniformly. There is more attraction at larger bag 
radii. The reason for this can be seen from Fig. 5, where the source function 
in Eq. ( [Ti"l) for R = 0.6 fm is plotted by the solid curve. dB/da and — g q a ip s 
are also shown by the dotted and the dashed curves, respectively. Due to 
dB/da the source function does not vanish even at large radii. Therefore, if 
the bag radius R is larger, <r(r) in Eq. becomes greater for large radii. 

This is shown in Fig. 6, where cr(r) is plotted for R — 0.6 and 0.8 fm. er(r) for 
R = 0.8 fm is greater than cr(r) for R = 0.6 fm at large radial region, which is 
the significant radial region. As a result, there is more attraction at large bag 
radii. 
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Fig. 5 : The source function is plotted by the solid curve. dB/dcr and — g&3p s are plotted by 
the dotted and the dashed curves, respectively. The coupling constants used are g% = gi = 2 
and = 3. 




r (fm) 

Fig. 6 : o"(r) is plotted with respect to the radial coordinate r by the solid and the dashed 
curves, respectively, for R = 0.6 and 0.8 fm. The coupling constants used are g q a = = 2 
and gf = 3. 

We made similar studies for other bag radii ranging from 0.6 to 1.0 fm and 
obtained similar results. Our results show that the change in the bag mass due 
to the coupling of the quarks with the mesons is not negligible. Such self-energy 
effects on the mass due to the meson coupling need to be taken into account in 
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choosing the parameters for the calculation of nuclear matter properties. The 
density-dependent bag constant can not only change the nucleon mass but also 
shift the bag radius. 

4 Summary 

We have applied the QMC model to a single nucleon. Recently the model has 
been often used to investigate explicit quark degrees of freedom in describing 
the nuclear matter. However, in the previous calculations the change in the 
mass of the bag due to the self-energy was ignored. Our calculations suggest 
that this change in the nucleon mass is not negligible, and thus the model 
parameters need to be modified to take this effect into consideration when 
they are used in nuclear matter calculations. 

5 Acknowledgement 

This work was supported by the Natural Sciences and Engineering Research 
Council of Canada. SWH acknowledges kind hospitality at TRIUMF and was 
partially supported by Sungkyun Faculty Research Fund and KOSEF Korea- 
Japan International Collaboration Program (986-0200-003-2). 

1. B. D. Scrot and J. D. Walecka, Adv. Nucl. Phys. 16, 1 (1986). 

2. P. A. M. Guichon, Phys. Lett. B 200, 235 (1988). 

3. A. Chodos, R. L. Jaffe, K. Johnson, C. B. Thorn, and V. Weisskopf, 
Phys. Rev. D 9, 3471 (1974); A. Chodos, R. L. Jaffe, K. Johnson, and 
C. B. Thorn, Phys. Rev. D 10, 2599 (1974). 

4. A. W. Thomas, Adv. Nucl. Phys. 13, 1 (1984). 

5. S. Fleck, W. Bentz, K. Shimizu, and K. Yazaki, Nucl. Phys. A510, 731 
(1990). 

6. K. Saito and A. W. Thomas, Phys. Lett. B 327 9 (1994). 

7. K. Saito and A. W. Thomas, Phys. Lett. B 335 17 (1994); 363 157 
(1995); Phys. Rev. C 51 2757 (1995); 52 2789 (1995). 

8. A. W. Thomas, A. Michels, A. W. Schreiber, and P. A. M. Guichon, 
Phys. Lett. B 233 43 (1989); K. Saito, A. Michels, and A. W. Thomas, 
Phys. Rev. C 46 R2149 (1992); K. Saito and A. W. Thomas, Nucl. 
Phys. A 574 659 (1994). 

9. H. Q. Song and R. K. Su, Phys. Lett. B 358 179 (1995). 

10. X. Jin and B. K. Jennings, Phys. Lett. B 374 13 (1996); Phys. Rev. C 
54 1427 (1996). 

11. J. D. Jackson, Classical Electrodynamics, 3rd., (John Wiley & Sons, Inc., 
New York, 1999) pp. 427. 



10 



